Process for estimating the channel in a OFDM communication system, and receiver for doing the same

ABSTRACT

A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots may include the steps of receiving, by a receiver, OFDM sample based subframes comprising a set of pilot signals, performing an iterative loop for the purpose of computing (M k , k k ) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(h k |y k , l k ) at time instant k, y k , and l k  respectively corresponding to the received signal and to a priori known information, the value of (M k , k k ) being computed as a function of (M k−1 , k k−1 ) as well as the values of the pilots received at instant k.

TECHNICAL FIELD

The invention relates to telecommunications and more particularly to a method for estimating the transmission channel vector in an Orthogonal Frequency Division Multiplexing (OFDM) communication system, and receiver for doing the same.

BACKGROUND ART

Orthogonal Frequency Division Multiplexing (OFDM) communication systems is a widely developed technique in the field of wireless digital communications thanks to the high possibilities offered by digital signal processing and Discrete Fourier Transform computations.

However, the efficiency of the process strongly depends on the estimation of the channel, which is achieved thanks to the communication of pilot signals within the subframe.

FIG. 1 illustrates the particular structure of the time-frequency OFDM grid, corresponding to a sub-frame, with the location of the pilots.

Many mechanisms are known in the art for computing an estimation of the channel vector h_(k) (k corresponding to the discrete instant k), including iterative algorithms deriving the particular knowledge of h_(k) from the knowledge of h_(k−1) and the new information brought be the pilots received at the instant k.

Generally speaking, the channel estimation is much more effective when one knows the law of variation of the channel—its time correlation—and when one may store the high volume of data derived from the received pilots so as to elaborate the most comprehensive knowledge of the channel estimation.

However, it is clear that the storage of a high volume of data increases the amount of storage required for the estimation process and moreover the amount of digital processing resources for processing such information.

For those reasons, many channel estimation processes are based on an iterative mechanism for computing the estimate ĥ_(k) of the vector channel only based on the preceding estimate ĥ_(k−1) and the newly received information. Such iterative mechanisms show clear limitations particularly when the channel shows a large delay spread prohibiting accurate interpolation between the OFDM subcarriers.

In such an environment, but not exclusively, there is a need for an effective iterative mechanism, which takes advantage of all the information being available at a given instant k−1 for computing the estimated channel ĥ_(k) at the instant k.

Such is the technical problem to be solved by the present invention

SUMMARY OF THE INVENTION

It is an object of the present invention to provide a process estimating a channel for an OFDM communication system which can be used in particular in a mobile handset.

It is a further object of the present invention to provide a process of channel estimation for an OFDM communication system which reduces the amount of data to be stored for computing an accurate estimate of the channel.

It is another object of the present invention to provide a channel estimation method, using the Maximum Entropy principle, which allows one to assign probability distributions for the channel, when the estimator only knows the last past inferred channel distribution and the new received pilot symbols.

Those and other objects are achieved by means of a process for estimating the channel in an OFDM communication system, with OFDM samples formed of a time-frequency grid comprising a set of pilots.

The process involves the steps of:

-   -   receiving OFDM subframes comprising a set of pilot signals;     -   performing an iterative loop for the purpose of computing         (M_(k), k_(k)) respectively corresponding to the covariance         matrix and the mean of the probability distribution of         p(h_(k)|y_(k), I_(k)) at time instant k, y_(k), and I_(k)         respectively corresponding to the received signal vector and to         a priori known information;         said value of (M_(k), k_(k)) being computed as a function of         (M_(k−1), k_(k−1)), λ being the time correlation (left to be         defined) between the true channels h_(k) and h_(k−1) as well as         the values of the pilots received at instant k.

More particularly, the process performs an estimation of the channel based on the mean k_(k) and variance M_(k) of the distribution of the channel estimation in accordance with the following formulation:

$\quad\left\{ \begin{matrix} {M_{k} = {{\lambda^{2}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k} = {{\lambda\; k_{k - 1}} + {\frac{1}{\sigma^{2}}M_{k}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}}} \right)}}}} \end{matrix} \right.$ With

-   N being the number of subcarriers in the OFDM symbol; -   k_(k) being the estimated channel at instant k; -   M_(k) being an N×N matrix representative of the variance of the     estimated channel at instant k; -   Q being a matrix initialized to a predetermined value; -   P_(k) is an N×N diagonal matrix comprising a “1” at the location     corresponding to a subcarrier at pilot position and a “0” elsewhere; -   I_(N) being the N×N identity matrix; -   σ² being the variance of the noise process; -   h′_(k) being equal to the received signal y_(k) divided pointwise by     the corresponding signal vector s_(k); -   λ corresponding to the channel time correlation between instant k−1     and k.

In one embodiment, the matrix Q is an a priori frequency covariance matrix for an L-tap uniform channel, i.e. Q=F_(L)F_(L) ^(H), with F_(L) the matrix formed by the L first columns of the N×N DFT matrix

Preferably, the time correlation factor (λ) corresponds to the correlation factor in the auto regressive model of the first order for the time series h_(k).

In one embodiment, the parameter M_(k) is forwarded to the base station as an information representative of the quality of the channel estimation.

Preferably, the process is used for the Long Term Evolution (LTE) standard.

A process is also provided, which can be used even without the knowledge of the time correlation variable λ.

In that case the process involves the following steps:

-   -   performing an iterative loop for the purpose of computing (M_(k)         ^((λ)), k_(k) ^((λ))) respectively corresponding to the variance         and the mean of the probability distribution of p(h_(k)|y_(k),         I_(k)) at time instant k, y_(k), and I_(k) respectively         corresponding to the received signal and to a priori known         information;         said value of (M_(k) ^((λ)), k_(k) ^((λ))) being computed as a         function of (M_(k−1) ^((λ)), k_(k−1) ^((λ))) for a predetermined         set of values of time correlation λ         {λ_(n)} with n=1 to M, for some positive integer M, between the         true channel h_(k) and h_(k−1), as well as the values of the         pilots received at instant k,

$\begin{matrix} {{{\hat{h}}_{k} = \frac{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}k_{k}^{(\lambda_{n})}}}{\frac{1}{N}{\sum\limits_{n = 1}^{N}\;{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}}}}}{With}\text{}{\quad\left\{ \begin{matrix} {M_{k}^{(\lambda)} = {{\lambda^{2}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k}^{(\lambda)} = {{\lambda\; k_{k - 1}^{(\lambda)}} + {\frac{1}{\sigma^{2}}M_{k}^{(\lambda)}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}^{(\lambda)}}} \right)}}}} \end{matrix} \right.}} & \; \\ {{And}\text{}\left\{ \begin{matrix} {{X(\lambda)} = \left( {I + {\frac{P_{k}}{\sigma^{2}}\left( {{\lambda^{2}M_{k - 1}^{(\lambda)}} + {\left( {1 - \lambda^{2}} \right)Q}} \right)}} \right)^{- 1}} \\ {{X(\lambda)} = {\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)^{H}{X(\lambda)}\frac{P_{k}}{\sigma^{2}}\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)}} \end{matrix} \right.} & \; \end{matrix}$ With k_(k) ^((λ)) being the estimated channel at instant k depending on λ,

-   M_(k) ^((λ)) being a matrix representative of the variance of the     estimated channel at instant k; depending on λ, -   Q being a matrix initialized to a predetermined value; -   P_(k) is a diagonal N×N matrix comprising a “1” at the location of a     pilot symbol; -   I_(N) being the N×N identity matrix; -   σ² being the variance of the noise; -   h′_(k) being equal to the received signal y_(k) pointwise divided by     the corresponding pilot s_(k).

At last, the invention encompasses a receiver for an OFDM communication system comprising:

-   -   a means for receiving a OFDM subframes comprising a set of         pilots signals;     -   a means for performing an iterative loop for the purpose of         computing (M_(k) , k_(k)) respectively corresponding to the         variance and the mean of the probability distribution of         p(h_(k)|y_(k), I_(k)) at time instant k, y_(k), and I_(k)         respectively corresponding to the received signal and to a         priori known information; said value of (M_(k), k_(k)) being         computed as a function of (M_(k−1), k_(k−1)), λ being the known         time correlation between the true channel estimate h_(k) and         h_(k−1), as well as the values of the pilots received at instant         k.

DESCRIPTION OF THE DRAWINGS

Other features of one or more embodiments of the invention will best be understood by reference to the following detailed description when read in conjunction with the accompanying drawings.

FIG. 1 illustrates a well known time-frequency grid of an OFDM subframe comprising pilot signals.

FIG. 2 illustrates a first embodiment of a process of channel estimation in accordance with the present invention.

FIG. 3 illustrates a second embodiment which does not require the knowledge of the channel time correlation factor.

FIGS. 4 and 5 illustrate simulation of the results obtained with the process.

DESCRIPTION OF THE PREFERRED EMBODIMENT

Some preliminary theoretical considerations (I) show to be appropriate before describing one embodiment of a channel estimation process (II)

I. Preliminary Theoretical Considerations

Channel estimation, along with most synchronization procedures, is an intricate information theoretic problem. Indeed, while Shannon (C. E. Shannon, “A mathematical theory of communications”, Bell System Technical Journal, vol. 27, no. 7, pp. 379-423, 1948) allows us to derive the capacity of a channel for which all synchronization parameters, plus the noise variance, are perfectly known, no such theory exists when the knowledge of some of these parameters is missing.

More precisely, for a scalar communication y=hx+n, x

, for some codebook

, if h is unknown, then the maximum a posteriori estimate for x is

$\begin{matrix} {\hat{x} = {\arg\;{\max\limits_{x \in X}{\int_{h}^{\;}{\int_{n}^{\;}{{p\left( {{x❘h},n} \right)}{p(h)}{p(n)}{\mathbb{d}h}{\mathbb{d}n}}}}}}} & (1) \end{matrix}$

which requires to have an a priori p(h) for h. But this a priori is too impractical to obtain and would require to know all possible channel realizations and their respective probability. As a consequence to this strong difficulty, most contributions in the synchronization field have provided various empirical models based on field observations in order either to give an expression to p(h) or, more practically, to propose good channel estimators ĥ to h.

The difficulty of handling estimation problems when little side information is available is treated by Jaynes, through the Bayesian probability field, thanks to the maximum entropy principle (MaxEnt) (E. T. Jaynes, “Information Theory and Statistical Mechanics”, Physical Review, APS, vol. 106, no. 4, pp. 620-630, 1957).

However MaxEnt does not allow to perform updates of probability when new information, such as new pilots in the channel estimation problem, is available. In this case, the complete set of past symbols along with a prior distribution for the channel h₀ at time t=0.

The inventors have discussed such question fro the OFDM framework in the article “A maximum entropy approach to OFDM channel estimation”, arxiv Preprint http://arxiv.org/pdf/0811.0778 when the channel delay spread, the channel time correlation and the signal-to noise ratio (SNR) are alternatively known or unknown. When these parameters are not perfectly known, MaxEnt provides channel estimates minimizing the estimate mean square error (MMSE estimates) that outperform classical estimates which use empirical (often erroneous) models. Recent contributions in the Bayesian probability field enable one to perform probability updates, in particular based minimum cross entropy considerations [See A. Caticha, “Lectures on Probability, Entropy and Statistical Physics”, arXiv:0808.0012v1 [physics.data-an], 2008].

1. OFDM System Model

Preliminary notations: in the following boldface lower case represent vectors, capital boldface characters denote matrices (I_(N) is the N×N identity matrix). The transposition operation is denoted (·)^(T). The Hermitian transpose is denoted (·)^(H). The operator diag(x) turns the vector x into a diagonal matrix. The symbol det(X) is the determinant of matrix X. The symbol E[·] denotes expectation. The Kronecker delta function is denoted δ_(x) that equals 1 if x=0 and equals 0 otherwise.

Consider a single cell OFDM system with N subcarriers. The cyclic prefix (CP) length is N_(CP) samples. In the timefrequency OFDM symbol grid, pilots are found in the symbol positions indexed by the function φ_(t)(n)

{0, 1} which equals 1 if a pilot symbol is present at subcarrier n, at symbol time index t, and 0 otherwise. The time-frequency grid is depicted in FIG. 1. Both data and pilots at time t are modeled by the frequency-domain vector S_(t)

C^(N) with pilot entries of amplitude |s_(t,k)|²=1. The transmission channel is denoted h_(t)

C^(N) in the frequency-domain with entries of variance E[|h_(t,k)|²]=1. The additive noise is denoted n_(t)

C^(N) with entries of variance E[|n_(t,k)|²]=σ². The noise process is modeled as Gaussian with independent and identically distributed (i.i.d.) entries, n˜

(0, σ²I_(N))

The time-domain representation of h_(t) is denoted v_(t) ε

^({dot over (L)}{dot over ( )})with L the channel length, i.e. the channel delay spread expressed in OFDM-sample unit. The frequency-domain received signal y_(t)

C^(N) is then: y _(t)=diag(h _(t))s _(t) +n _(t)   (2)

We will also denote, ∀k ε {1, . . . , N} h′ _(k) =y _(k) /s _(k) =h _(k) +n _(k) /s _(k)

and h′=(h′₁, . . . , h′_(N))^(T) (here, the time index t is implicit).

The channel h_(t) evolves in time with coherence time function λ(τ) such that, independently of the channel delay spread index

$\begin{matrix} {{E\left\lbrack {v_{i,t}v_{i,{t + \tau}}^{*}} \right\rbrack} = \frac{\lambda(\tau)}{L}} & (3) \end{matrix}$

This model is justified by maximum entropy considerations (see e.g. the reference R. Couillet, M. Debbah, “A maximum entropy approach to OFDM channel estimation”, arxiv Preprint http://arxiv.orq/pdf/0811.0778) and leads for p(h_(t)|h_(t+τ)) to be assigned a Gaussian distribution of mean λ(τ)h_(t+τ) and variance

$\begin{matrix} {{{E\left\lbrack {h_{t}h_{t + \tau}^{H}} \right\rbrack} = {\left( {1 - {\lambda(\tau)}^{2}} \right)Q}}{With}} & (4) \\ {Q_{nm} = {{E\left\lbrack {\sum\limits_{k = 0}^{L - 1}{\sum\limits_{l = 0}^{L - 1}{v_{k}v_{l}^{*}{\mathbb{e}}^{{- 2}{\pi\mathbb{i}}\frac{{kn} - {lm}}{N}}}}} \right\rbrack} = {\frac{1}{L}{\sum\limits_{k = 0}^{L - 1}{\mathbb{e}}^{{- 2}{\pi\mathbb{i}}\; k\frac{n - m}{N}}}}}} & (5) \end{matrix}$

This work aims at providing MMSE estimates ĥ_(t)

C^(N) of the vector ht when the receiver, at time t is only aware of:

(i) the last stream of pilots received at some time t_(p)≦t,

(ii) the last density function p(h_(t0)) assigned to h_(t0). All other system parameters, e.g. channel delay spread, channel coherence time, SNR, are either known or unknown. This will be clearly stated in the following and will be denoted I.

2. ME Principle and Application to Channel Estimation.

The fundamental difference between orthodox and Bayesian probabilities lies in the existence of prior distributions in the Bayesian realm. In the present situation, when estimating h_(t) from the observation y_(t), orthodox probabilists decide for a maximum likelihood estimator ĥ_(t)=arg max_(ht) p(h_(t)|h′_(t)), following a given model linking h′_(t) to h_(t). Bayesian probabilists would note that p(h _(t) |h′ _(t) , I)=p(h′ _(t) |h _(t) , I)p(h _(t) |I)/p(h′ _(t) |I). p(h _(t) |h′ _(t) , I)=p(h′ _(t) |h _(t) , I)p(h _(t) |I)/p(h′ _(t) |I),

where I is the side information of the estimator. The probability p(h_(t)|I) is the subjective prior distribution for h_(t), which is immaterial in orthodox probability.

To derive p(h_(t)|I) based on prior statistical knowledge on h_(t), e.g. first order moments, Jaynes introduces the maximum entropy principle in the reference E. T. Jaynes, “Information Theory and Statistical Mechanics”, Physical Review, APS, vol. 106, no. 4, pp. 620-630, 1957.

The latter allows to assign a distribution to p(h_(t)|I) which satisfies the statistical constraints given by I and which minimizes the risk of making unwanted assumptions, by maximizing the entropy of p(h_(t)|I).

When it comes to update probability assignments, Caticha proposes an extension of the maximum entropy principle, namely the minimum cross entropy principle (ME) in A. Caticha, “Lectures on Probability, Entropy and Statistical Physics”, arXiv:0808.0012v1 [physics.data-an], 2008.

When p(h_(t) |I ₁) has been assigned for some side information I₁, and new cogent information I₂ is later available, then the ME principle consists in assigning to p(h_(t)|I₂) the distribution

$\begin{matrix} {{{p\left( {h_{t}❘I_{2}} \right)} = {\arg\;{\min\limits_{q}\;{S\left\lbrack {q,{p\left( {h_{t}❘I_{1}} \right)}} \right\rbrack}}}}{Where}} & (6) \\ {{S\left\lbrack {q,p} \right\rbrack} = {\int_{\;}^{\;}{{q(x)}{\log\left( \frac{p(x)}{q(x)} \right)}{\mathbb{d}x}}}} & (7) \end{matrix}$

It can therefore be derived a method based on a minimal update requirement, which in essence assigns to p(h_(t)|I₂) the unique distribution which minimizes the changes brought to p(h_(t)|I₁) while satisfying the new constraints given by I₂.

In the following, additional side information on h_(t) (which possibly varies over time) will come from new available pilots at later time positions.

3. dynamic Optimal Channel Estimation for OFDM with Perfect System Parameter Knowledge

Channel estimation is performed at different time instants t_(k), k ε

.

Denote I_(k) the knowledge at time t_(k). Since memory restrictions impose to discard past, we decide here only to consider at time t_(k) the last received pilot data symbols, the last assigned probability p(h_(k)|I_(k−1)), h_(k)=h_(tk) and the supposedly known time correlation λ between time t and the past channel h_(k−1). Assume prior assigned distribution p(h_(k−1)|I_(k−1)) at time index k−1. We have in general

$\begin{matrix} {p\left( {{h_{k}❘y_{k}},I_{k}} \right)} & (8) \\ {= \frac{{p\left( {{y_{k}❘h_{k}},I_{k}} \right)} \cdot {p\left( {h_{k}❘I_{k}} \right)}}{p\left( {y_{k}❘I_{k}} \right)}} & (9) \\ {= \frac{{p\left( {{y_{k}❘h_{k}},I_{k}} \right)} \cdot {\int{{{p\left( {{h_{k}❘h_{k - 1}},I_{k}} \right)} \cdot {p\left( {h_{k - 1}❘I_{k}} \right)}}{\mathbb{d}h_{k - 1}}}}}{p\left( {y_{k}❘I_{k}} \right)}} & (10) \end{matrix}$

Using Caticha's ME principle and set p(h_(k−1)|I_(k)) to the previous p(h_(k−1)|y_(k−1), I_(k−1)). The reason lies in the minimal update principle: if no additional information is given in I_(k), compared to I_(k−1), then p(h_(k−1)|y_(k−1), I_(k−1)) is the distribution q that minimizes the cross-entropy S[q, p(h_(k−1)|y_(k−1), I_(k−1))]²

Let us now perform a recursive reasoning over the channel estimates at time indexes k ε

. Assume that p(h_(k−1)|y_(k−1), I_(k−1)) is Gaussianz,

(k_(k−1), M_(k−1)).

One may show that this implies p(h_(k)|y_(k), I_(k)) is still Gaussian. This will therefore be denoted

(k_(k), M_(m) ⁻¹)

We have

$\begin{matrix} {\mspace{79mu}{p\left( {{h_{k}❘y_{k}},I_{k}} \right)}} & (11) \\ {\mspace{79mu}{= {\alpha_{1}{{p\left( {{y_{k}❘h_{k}},I_{k}} \right)} \cdot {\int{{{p\left( {{h_{k}❘h_{k - 1}},I_{k}} \right)} \cdot {p\left( {h_{k - 1},I_{k}} \right)}}{\mathbb{d}h_{k - 1}}}}}}}} & (12) \\ {= {\lim\limits_{\overset{\sim}{Q}\longrightarrow Q}{{\mathbb{e}}^{{({h_{k} - h_{k}^{\prime}})}^{H}\frac{P_{k}}{\sigma^{2}}{({h_{k} - h_{k}^{\prime}})}}{\int{{\mathbb{e}}^{{({h_{k} - {\lambda\; h_{k - 1}}})}^{H}\frac{{\overset{\sim}{Q}}^{- 1}}{1 - \lambda^{2}}{({h_{k} - {\lambda\; h_{k - 1}}})}} \times \alpha_{2}{\mathbb{e}}^{{({h_{k - 1} - k_{k - 1}})}^{H}{M_{k - 1}^{- 1}{({h_{k - 1} - k_{k - 1}})}}}{\mathbb{d}h_{k - 1}}}}}}} & (13) \end{matrix}$

where the {tilde over (Q)}'s are taken from a set of invertible matrices in the neighborhood of Q, and the α_(i)'s are constants.

First we need to write the exponents of the Gaussian products in the integrand in a single Gaussian exponent form of the vector h_(k−1) times a constant independent of h_(k−1). By expansion and simplification, this is

$\begin{matrix} {{{{\left( {h_{k} - {\lambda\; h_{k - 1}}} \right)^{H}\frac{{\overset{\sim}{Q}}^{- 1}}{1 - \lambda^{2}}\left( {h_{k} - {\lambda\; h_{k - 1}}} \right)} + {\left( {h_{k - 1} - k_{k - 1}} \right)^{H}{M_{k - 1}^{- 1}\left( {h_{k - 1} - k_{k - 1}} \right)}}} = {{\left( {h_{k - 1} - 1} \right)^{H}{N\left( {h_{k - 1} - 1} \right)}} + {C\left( h_{k} \right)}}}{With}\left\{ \begin{matrix} {N = {\frac{\lambda^{2}{\overset{\_}{Q}}^{- 1}}{1 - \lambda^{2}} + M_{k - 1}^{- 1}}} \\ {1 = {N^{- 1}\left( {{\frac{\lambda}{1 - \lambda^{2}}{\overset{\sim}{Q}}^{- 1}h_{k}} + {M_{k - 1}^{- 1}k_{k - 1}}} \right)}} \\ {C = {{h_{k}^{H}\frac{{\overset{\_}{Q}}^{- 1}}{1 - \lambda^{2}}h_{k}} + {k_{k - 1}^{H}M_{k - 1}^{- 1}k_{k - 1}} - {\left( {h_{k} + {\frac{1 - \lambda^{2}}{\lambda}\overset{\sim}{Q}M_{k - 1}^{- 1}}} \right)^{H}\frac{\lambda^{2}}{\left( {1 - \lambda^{2}} \right)^{2}}{{\overset{\sim}{Q}}^{- 1}\left( {\frac{\lambda^{2}{\overset{\_}{Q}}^{- 1}}{1 - \lambda^{2}} + M_{k - 1}^{- 1}} \right)}^{- 1}{{\overset{\sim}{Q}}^{- 1}\left( {h_{k} + {\frac{1 - \lambda^{2}}{\lambda}\overset{\sim}{Q}M_{k - 1}^{- 1}k_{k - 1}}} \right)}}}} \end{matrix} \right.} & (14) \end{matrix}$

The integral (13) is then a constant times e^(C), which depends on h_(k). The term C must then be written again into a quadratic expression of h_(k). This is

$\begin{matrix} {{C = {{\left( {h_{k} - j} \right)^{H}{R\left( {h_{k} - j} \right)}} + B}}{with}} & (15) \\ \left\{ \begin{matrix} {R = \left( {{\lambda^{2}M_{k - 1}} + {\left( {1 - \lambda^{2}} \right)\overset{\sim}{Q}}} \right)^{- 1}} \\ {j = {\lambda\; k_{k - 1}}} \\ {B = 0} \end{matrix} \right. & (16) \end{matrix}$

Together with the term outside the integral (13), this is p(h _(k) |y _(k) , I _(k))=α·e ^((h) ^(k) ^(−k) ^(k) ⁾ ^(H) ^(M) ^(k) ⁻¹ ^((h) ^(k) ^(−k) ^(k) ⁾   (17)

With α=(∫p(h_(k)|y_(k), I_(k))dh_(k))⁻¹.

Finally, after some arithmetic derivation, in the limit {tilde over (Q)}→Q,

$\begin{matrix} {\quad{\quad\left\{ \begin{matrix} {M_{k} = {{\lambda^{2}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k} = {{\lambda\; k_{k - 1}} + {\frac{1}{\sigma^{2}}M_{k}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}}} \right)}}}} \end{matrix} \right.}} & (18) \end{matrix}$

And the MMSE estimator ĥ_(k) for the channel at time index k is the first order moment of a Gaussian distribution centered in k_(k) which is ĥ_(k)=k_(k). At initial time instant t₀, if nothing but the channel delay spread L is known, M₀=Q from MaxEnt, and k₀=0. Therefore, we prove by the above recursion that, under this state of initial knowledge, for all k ε

p(h_(k)|I_(k)) is Gaussian with mean k_(k) and variance M_(k), and ĥ_(k)=k_(k). Note that, while regularized inverses of Q were used along the derivations, the final formulas are properly conditioned with respect to Q

The process which will be described hereafter will take advantage of this particular observation, and will achieve an iterative mechanism of channel estimation, based on an iterative update of both the mean k_(k) and variance M_(k) of the distribution of the channel estimation.

Before entering into the particular detail of the process, the inventors have further discovered that the channel estimation can even be achieved without perfect system parameter knowledge 4. Dynamic Optimal Channel Estimation for OFDM with Perfect System Parameter Knowledge

In practical applications, contrary to what was stated above, the different parameters λ, σ² and L especially, might not be perfectly known. For simplicity we assume those parameters are constant over the duration of the channel estimation process. Following the MaxEnt principle, these parameters must be assigned an a priori distribution. Let us focus on the time correlation λ, which is typically the most difficult parameter to track. In this respect, one has p(h _(k) |y _(k) , I _(k))=∫p(h _(k) |y _(k) , λI _(k))p(λ|y _(k) , I _(k))dλ  (19)

Since y_(k) cannot bring alone any cogent information on λ, p(λ|y_(k), I_(k))=p(λ|I_(k)). The probability p(h_(k)|y_(k), λ, I_(k)) was computed above and is given by the right-hand side of Equation (17), in which a depends on A and must therefore be made explicit.

Further computation leads to p(h _(k) |y _(k) , I _(k))=∫p(λ|I _(k))α(λ)e ^((h) ^(k) ^(−k) ^(k) ^((λ)) ⁾ ^(H) ^((M) ^(k) ^((λ)) ⁾ ⁻¹ ^((h) ^(k) ^(−k) ^(k) ^((λ)) ⁾ dλ  (20)

with M_(k) ^((λ)) and k_(k) ^((λ)) given by Equation (18) above for the λ in question and

$\begin{matrix} {{{\alpha(\lambda)} = {{\beta\mathbb{e}}^{- {x{(\lambda)}}}{\det\left\lbrack {X(\lambda)} \right\rbrack}}}{With}} & (21) \\ \left\{ {{{\begin{matrix} {{X(\lambda)} = \left( {I + {\frac{P_{k}}{\sigma^{2}}\left( {{\lambda^{2}M_{k - 1}^{(\lambda)}} + {\left( {1 - \lambda^{2}} \right)Q}} \right)}} \right)^{- 1}} \\ {{x(\lambda)} = {\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)^{H}{X(\lambda)}\frac{P_{k}}{\sigma^{2}}\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)}} \end{matrix}{And}\beta} = \left( {\int{{p\left( {{h_{k}❘y_{k}},I_{k}} \right)}{\underset{\bigwedge}{\mathbb{d}}h_{k}}}} \right)^{- 1}},{{independent}\mspace{14mu}{of}\mspace{14mu}{\lambda.}}} \right. & (22) \end{matrix}$

The MMSE estimate ĥ_(k) is then the weighted sum:

$\begin{matrix} {{\hat{h}}_{k} = \frac{\int{{p\left( {\lambda ❘I_{k}} \right)}{\mathbb{e}}^{- {x{(\lambda)}}}{\det\left\lbrack {X(\lambda)} \right\rbrack}k_{k}^{(\lambda)}{\mathbb{d}\lambda}}}{\int{{p\left( {\lambda ❘I_{k}} \right)}{\mathbb{e}}^{- {x{(\lambda)}}}{\det\left\lbrack {X(\lambda)} \right\rbrack}{\mathbb{d}\lambda}}}} & (23) \end{matrix}$

This integral is however very involved. In practice, it must be broken into a finite sum over a set of potential values for λ. Denoting S this set and |S| its cardinality, the recursive algorithm that provides the successive estimates ĥn_(k) k=1, . . . , K, requires that at every step, the values for M_(k) ^((λ)) and k_(k) ^((λ)), λ

S are kept in memory.

One may then describe the preferred embodiment of a channel estimation process which is based on the underlying above described principle.

II. Preferred Embodiment of the Channel Estimation Process for a OFDM Communication System

It has been demonstrated how the knowledge of the variance and estimate M_(k) ^((λ)) and k_(k) ^((λ)) at k−1 can be advantageously used for achieving an iterative estimation method based on the update mechanism shown in formulas (18). With the hypothesis of a Gaussian channel, the information collected by the estimate at instant k−1 and the variance shows to be minimal and, therefore, allows to reduce the quantity of data to be stored within the memory since the variance will only need to be represented by an N×N matrix.

Furthermore, while the processes described below will be discussed in the frequency domain, where a matrix shows to be an N×N matrix which might have a significant size, the computations below may also be performed in the time domain what will result in a significant reduction of the size of the matrix.

With respect to FIG. 2, there is now described one embodiment of a process which assumes the knowledge of the value of λ

In a step 21, the process starts with an initialization step. At time t=0, the values of M₀ and k₀ are respectively set to Q and 0.

Then, in a step 22, the process enters into a loop;

Then, in a step 23, the process receives an OFDM symbol and may extract the received signal and the pilots therein contained.

Then, in a step 24, the process computes an estimate value of the variance of the probability distribution of p(h_(k)|y_(k), I_(k)) at time instant k, y_(k), and I_(k) respectively corresponding to the received signal and to a priori known information, in accordance with the following formula:

$M_{k} = {{\lambda^{2}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}$ With

-   N being the number of sub carriers in the OFDM symbol; -   k_(k) being the estimated channel at instant k; -   M_(k) being a N×N matrix representative of the variance of the     estimated channel at instant k; -   Q being a matrix initialized to a predetermined value; -   P_(k) is an N×N diagonal matrix comprising a “1” at the location of     a pilot symbol and a “0” elsewhere; -   I_(N) being the N×N identity matrix; -   λ corresponding to the known channel time correlation between     instant k−1 and k.

Then, in a step 25the process computes an estimate value of the channel k_(k) in accordance with the formula:

$k_{k} = {{\lambda\; k_{k - 1}} + {\frac{1}{\sigma^{2}}M_{k}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}}} \right)}}}$ With

-   σ² being the variance of the noise; -   h′_(k) being equal to the received signal y_(k) pointwise divided by     the corresponding pilot s_(k);

Then, in a step 26, the process issues the estimate value ĥ_(k)=k_(k)

Then, in a step 27, the process loops back to step 22 in order to process the next OFDM symbol.

With respect to FIG. 3, there is illustrated a second embodiment of a process of channel estimation which takes no assumption on the value of λ.

In a step 31, the process starts with an initialization step similar to step 21 above.

In a step 32 the process enters a loop.

In a step 33, the process receives an OFDM sample.

Then, in step 34, the process enters a new loop for processing a predetermined set of possible coefficients number of values (λ₁ λ₂ λ₃ . . . λ_(n) . . . λ_(M)) with n=1 to M. To each value λ_(n) one associates a couple (M_(k) ^((λ)), k_(k) ^((λ)), which will be calculated in accordance with the same formula as above.

Indeed, in steps 35 and 36 the process successively computes M_(k) ^((λ)) and k_(k) ^((λ)) in accordance with the formulas:

$\quad\left\{ \begin{matrix} {M_{k}^{(\lambda)} = {{\lambda^{2}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k}^{(\lambda)} = {{\lambda\; k_{k - 1}^{(\lambda)}} + {\frac{1}{\sigma^{2}}M_{k}^{(\lambda)}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}^{(\lambda)}}} \right)}}}} \end{matrix} \right.$ With

-   N being the number of sub carriers in the OFDM symbol; -   k_(k) being the estimated channel at instant k; -   M_(k) being a N×N matrix representative of the variance of the     estimated channel at instant k; -   Q being a matrix initialized to a predetermined value; -   P_(k) is an N×N diagonal matrix comprising a “1” at the location of     a pilot symbol and a “0” elsewhere; -   I_(N) being the N×N identity matrix; -   λ corresponding to the known channel time correlation between     instant k−1 and k.

Then, in a step 37, a test is considered in order to determine whether a subsequent value λ_(n+1) is to be processed.

If so, the process returns to step 34 above for the purpose of computing that new value and, conversely, if all the values of (λ₁ λ₂ λ₃ . . . λ_(n) . . . λ_(M)) have been processed, then the process proceeds to step 38 where it computes the estimate ĥ_(k) of the channel from a weighted sum of the latter, as follows:

${\hat{h}}_{k} = \frac{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}k_{k}^{(\lambda_{n})}}}{\frac{1}{N}{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}}}}$ With $\left\{ \begin{matrix} {{X(\lambda)} = \left( {I + {\frac{P_{k}}{\sigma^{2}}\left( {{\lambda^{2}M_{k - 1}^{(\lambda)}} + {\left( {1 - \lambda^{2}} \right)Q}} \right)}} \right)^{- 1}} \\ {{x(\lambda)} = {\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)^{H}{X(\lambda)}\frac{P_{k}}{\sigma^{2}}\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)}} \end{matrix} \right.$ With k_(k) ^((λ)) being the estimated channel at instant k depending on λ,

M_(k) ^((λ)) being a matrix representative of the variance of the estimated channel at instant k assuming λ time correlation.

The process then proceeds to step 39 which loops back to step 32 for the purpose of processing a new OFDM symbol.

Simulations and Results

In this section, we provide simulation plots to compare, at time t, the minimal channel estimation update method against (i) the one-dimensional MMSE, taken at time t, which takes only into account the last past pilot symbols and uses a fixed empirical covariance matrix, (ii) an typical optimum two-dimensional MMSE, with K=4 pilot time indexes. The OFDM DFT size is N=64, the channel length L=6 is known to the receiver, the vehicular speed is v=120 km/h, pilot sequences are transmitted every 0.29 ms (as in 3GPPLTE [7]), and the induced Jake's time correlation λ between t and the past pilot sequence arrival time is known to the receiver. In scenario (ii), all K past received pilot sequences and time correlations are perfectly known. The channel time correlation model is a K-order autoregressive model following the reference K. E. Baddour and, N. C. Beaulieu, “Autoregressive modeling for fading channel simulation,” IEEE Transaction on Wireless Communications, July 2005.

A performance comparison is proposed in FIG. 5. We notice here that the minimal update algorithm does not show significant performance decay compared to the optimal two-dimensional MMSE estimator, while the one-dimensional MMSE estimator, also relying on the last past pilot sequence, shows large performance impairment.

In FIG. 3, with the same assumptions as previously, we consider the hypothesis where the vehicular speed v is a priori known to be (with equal probabilities) either 5, 50, 120 km/h. The performance is compared against the optimal 2-D algorithm where v is known but erroneously estimated (v=5, 50, 100 km/h). It is observed that, again, even when λ, or equivalently v, is a priori unknown, the Bayesian minimal update framework manages to ideally recover the channel with no performance decay. On the opposite, when λ is erroneously estimated, the performance decay of the optimal estimator might be dramatic. 

The invention claimed is:
 1. A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots, the process comprising the steps of: receiving, by a receiver, OFDM sample based subframes comprising a set of pilot signals; performing an iterative loop for the purpose of computing (M_(k), k_(k)) respectively corresponding to the covariance matrix and the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k, y_(k), and l_(k) respectively corresponding to the received signal and to a priori known information and h_(k) being a true channel at time instant k; and said value of (M_(k), k_(k)) being computed as a function of (M_(k−1), k_(k−1)) as well as the values of the pilots received at instant k.
 2. The process for estimating the channel in an OFDM communication system according to claim 1 further comprising the steps of: receiving OFDM sample based subframes comprising a set of pilots signals; performing an estimation of the channel based on the mean k_(k) and variance M_(k) of the distribution of the true channel h_(k) in accordance with the following formulation: $\quad\left\{ \begin{matrix} {M_{k} = {{\lambda^{2}\left( {M_{k - 1} + {\frac{1 - \lambda^{3}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k} = {{\lambda\; k_{k - 1}} + {\frac{1}{\sigma^{2}}M_{k}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}}} \right)}}}} \end{matrix} \right.$ with N being the number of sub carriers in the OFDM symbol; k_(k) being the mean of the inferred probability distribution p(h_(k)|y_(k), l_(k)) at time instant k; M_(k) being a N×N matrix representative of the variance of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k; Q being a matrix initialized to a predetermined value; P_(k) is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere; I_(N) being the N×N identity matrix; λ corresponding to the known channel time correlation between instant k−1 and k; σ² being the variance of the noise; and h′_(k) being equal to the received signal y_(k) pointwise divided by the corresponding pilot s_(k).
 3. The process according to claim 2 wherein Q is an a priori frequency covariance matrix for an L-tap uniform channel.
 4. The process according to claim 2 wherein the time correlation factor (λ) corresponds to the correlation factor in the auto regressive model of the first order of h_(k).
 5. The process according to claim 1 wherein said parameter M_(k) is forwarded to the base station as an information representative of the quality of the channel estimation.
 6. The process according to claim 1 adapted to the 3GPP Long Term Evolution (LTE) standard.
 7. A process for estimating a channel in an OFDM communication system, wherein an OFDM sample is based on a time-frequency grid comprising a set of pilots, the process comprising the steps of: receiving, by a receiver, OFDM sample based subframes comprising a set of pilots signals; performing an iterative loop for the purpose of computing (M_(k) ^((λ)), k_(k) ^((λ))) respectively corresponding to the variance and the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k, y_(k), and l_(k) respectively corresponding to the received signal and to a priori known information; said value of (M_(k) ^((λ)), k_(k) ^((λ))) being computed as a function of (M_(k−1) ^((λ)), k_(k−1) ^((λ))) for λspanning a predetermined set of values {λ_(n)} with n=1 to M between true channel h_(k) and h_(k−1) respectively at time instant k and k−1, as well as the values of the pilots received at instant k, $\begin{matrix} {{{\hat{h}}_{k} = \frac{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}k_{k}^{(\lambda_{n})}}}{\frac{1}{N}{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}}}}}{With}\left\{ \begin{matrix} {M_{k}^{(\lambda)} = {{\lambda^{2}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k}^{(\lambda)} = {{\lambda\; k_{k - 1}^{(\lambda)}} + {\frac{1}{\sigma^{2}}M_{k}^{(\lambda)}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}^{(\lambda)}}} \right)}}}} \end{matrix} \right.} & \; \\ {{And}\left\{ \begin{matrix} {{X(\lambda)} = \left( {I + {\frac{P_{k}}{\sigma^{2}}\left( {{\lambda^{2}M_{k - 1}^{(\lambda)}} + {\left( {1 - \lambda^{2}} \right)Q}} \right)}} \right)^{- 1}} \\ {{x(\lambda)} = {\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)^{H}{X(\lambda)}\frac{P_{k}}{\sigma^{2}}\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)}} \end{matrix} \right.} & \; \end{matrix}$ with k_(k) ^((λ)) being the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at instant k depending on λ, . . . M_(k) ^((λ)) being a matrix representative of the variance of the estimated channel at instant k; depending on λ, Q being a matrix initialized to a predetermined value; P_(k) is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol; I_(N) being the N×N identity matrix; σ² being the variance of the noise; h′_(k) being equal to the received signal y_(k) pointwise divided by the corresponding pilot s_(k;) λ corresponding to the known channel time correlation between instant k−1 and k; H is the Hermitian transpose; and det[X] is the determinant of the matrix X.
 8. The process according to claim 7 wherein Q is an a priori frequency correlation matrix for an L-tap uniform channel.
 9. The process according to claim 7 wherein the set of time correlation factors {λ_(n)} corresponds to hypothetical values for the time correlation parameter λ.
 10. The process according to claim 7 wherein said parameter M_(k) is forwarded to the base station as an information representative of the quality of the channel estimation.
 11. The process according to claim 7 wherein it is adapted to the Long Term Evolution (LTE) standard.
 12. A receiver for an OFDM communication system comprising: a means for receiving a OFDM sample based subframes comprising a set of pilots signals; a means for performing an iterative loop for the purpose of computing (M_(k) , k_(k)) respectively corresponding to the variance and the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k, y_(k), and l_(k) respectively corresponding to the received signal and to a priori known information; and said value of (M_(k), k_(k)) being computed as a function of (M_(k−1), k_(k−1)) as well as the values of the pilots received at instant k, with h_(k) corresponding to the channel at the kth sampled time instant.
 13. A receiver according to claim 12, said receiver being incorporated to a mobile communication device for an LTE communication system.
 14. The receiver according to claim 12, further comprising: a means for performing an estimation of the channel based on the mean k_(k) and variance M_(k) of the distribution of the channel estimation in accordance with the following formulation: $\quad\left\{ \begin{matrix} {M_{k} = {{\lambda^{2}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k} = {{\lambda\; k_{k - 1}} + {\frac{1}{\sigma^{2}}M_{k}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}}} \right)}}}} \end{matrix} \right.$ with N being the number of sub carriers in the OFDM symbol; k_(k) being the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k; M_(k) being a N×N matrix representative of the variance of the probability distribution of p(h_(k)|y_(k), l_(k)) at instant k; Q being a matrix initialized to a predetermined value; P_(k) is an N×N diagonal matrix comprising a “1” at the location of a pilot symbol and a “0”elsewhere; I_(N) being the N×N identity matrix; and λ corresponding to the known channel time correlation between instant k−1 and k; σ² being the variance of the noise; and h′_(k) being equal to the received signal y_(k) pointwise divided by the corresponding pilot s_(k).
 15. A receiver for an OFDM communication system comprising a means for estimating the channel and characterized in that it further comprises: a means for receiving a OFDM sample based subframes comprising a set of pilots signals; a means for performing an iterative loop for the purpose of computing (M_(k) ^((λ)), k_(k) ^((λ))) respectively corresponding to the variance and the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at time instant k, y_(k), and l_(k) respectively corresponding to the received signal and to a priori known information; said value of (M_(k) ^((λ)), k_(k) ^((λ))) being computed as a function of (M_(k−1), k_(k−1)) for a predetermined set of values of time correlation (λ_(n)) with n=1 to M between the true channels h_(k) and h_(k−1), with h_(k) corresponding to the channel at the kth sampled time instant,as well as the values of the pilots received at instant k, $\begin{matrix} {{{\hat{h}}_{k} = \frac{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}k_{k}^{(\lambda_{n})}}}{\frac{1}{N}{\sum\limits_{n = 1}^{N}{{\mathbb{e}}^{- {x{(\lambda_{n})}}}{\det\left\lbrack {X\left( \lambda_{n} \right)} \right\rbrack}}}}}{with}\left\{ \begin{matrix} {M_{k}^{(\lambda)} = {{\lambda^{2}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)} \times \left( {{\frac{\lambda^{2}}{\sigma^{2}}{P_{k}\left( {M_{k - 1}^{(\lambda)} + {\frac{1 - \lambda^{2}}{\lambda^{2}}Q}} \right)}} + I_{N}} \right)^{- 1}}} \\ {k_{k}^{(\lambda)} = {{\lambda\; k_{k - 1}^{(\lambda)}} + {\frac{1}{\sigma^{2}}M_{k}^{(\lambda)}{P_{k}\left( {h_{k}^{\prime} - {\lambda\; k_{k - 1}^{(\lambda)}}} \right)}}}} \end{matrix} \right.} & \; \\ {{and}\left\{ \begin{matrix} {{X(\lambda)} = \left( {I + {\frac{P_{k}}{\sigma^{2}}\left( {{\lambda^{2}M_{k - 1}^{(\lambda)}} + {\left( {1 - \lambda^{2}} \right)Q}} \right)}} \right)^{- 1}} \\ {{x(\lambda)} = {\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)^{H}{X(\lambda)}\frac{P_{k}}{\sigma^{2}}\left( {{\lambda\; k_{k - 1}^{(\lambda)}} - h_{k}^{\prime}} \right)}} \end{matrix} \right.} & \; \end{matrix}$ With k_(k) ^((λ)) being the mean of the probability distribution of p(h_(k)|y_(k), l_(k)) at instant k depending on λ, M_(k) ^((λ)) being a matrix representative of the variance of the estimated channel at instant k; depending on λ, Q being a matrix initialized to a predetermined value; P_(k) is a diagonal N×N matrix comprising a “1” at the location of a pilot symbol; I_(N) being the N×N identity matrix; σ² being the variance of the noise; and h′_(k) being equal to the received signal y_(k) pointwise divided by the corresponding pilot s_(k;) λ corresponding to the known channel time correlation between instant k−1 and k; H is the Hermitian transpose; and det[X] is the determinant of the matrix X. 